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We present a theoretical study of the quantum magnetocapacitance of an ultrathin topological insulator film 
in an external magnetic field. The study is undertaken to investigate the interplay of the Zeeman interaction 
with the hybridization between the upper and lower surfaces of the thin film. Determining the density of 
states, we find that the electron-hole symmetry is broken when the Zeeman and hybridization energies are 
varied relative to each other. This leads to a change in the character of the magnetocapacitance at the charge 
neutrality point. We further show that in the presence of both Zeeman interaction and hybridization the 
magnetocapacitance exhibits beating at low and splitting of the Shubnikov de Haas oscillations at high 
perpendicular magnetic field. In addition, we address the crossover from perpendicular to parallel magnetic 
field and find consistency with recent experimental data. 

Topological insulators (TIs) form a new class of compounds with strong spin-orbit interaction 18 . These 
materials possess an energy gap in the bulk (like insulators) and conducting electronic states at the surface 
(like metals). The nature of the surface states depends on the dimensions of the TP" 16 : Two-dimensional 
surface states have been revealed in transport experiments for three-dimensional TIs. The charge carriers in such 
surface states form a two-dimensional gas of Dirac fermions with a single Dirac cone energy spectrum. Dirac 
fermions have been observed in angle resolved photo-emission experiments and quantization of their Landau 
levels (LLs) has been confirmed by scanning tunneling spectroscopy 17 21 . TI thin films provide a new degree of 
freedom which is the thickness of the film. Effects of the crossover from 3D to 2D topological surface states as well 
as hybridization of the states at the two surfaces of the thin film have been studied for varying film thicknesses and 
found to become important when the thickness is less than about 5 nm M . The Dirac fermions in TIs share 
several features with those of graphene 27 ' 28 . However, there are important differences between graphene and the 
surface states of a TI. Unlike graphene, there is an odd number of not spin degenerate Dirac points at any surface. 
In addition, TIs exhibit a unique spin-momentum locking in the surface states. 

Observation of magnetic quantum oscillations, Shubnikov de Haas (SdH) and de Haas van Alphen 29,30 , in the 
transport experiments has resulted in theoretical work on magnetization 29 and Berry phase 3031 features of TI 
surface Dirac fermions. Furthermore, much attention is paid to magnetically doped TIs due to their exceptional 
properties and potential spintronics applications 10,24 . In this context, capacitance measurements are among the 
most important tools for studying the electronic properties. They can be effectively used to probe the thermodyn- 
amic density of states (DOS). Although the focus in TI research has been on transport properties, insight into the 
fundamental electronic properties and device physics still calls for knowledge about the capacitance-voltage (C- 
V) characteristics of the system. Results have been reported for carbon nanotubes, graphene nanoribbons, and 
mono- and bilayer graphene systems 32 37 . In view of this, attention is now being paid to electrostatic properties, 
such as the magnetocapacitance of TI thin film devices. Furthermore, for improve the performance of field effect 
transistors 38 , the potential of TI thin films as channel materials is creating much excitement. This is due to the 
excellent intrinsic transport features as well as the possibility of patterning device structures within top-down 
lithographical techniques. The present work aims at determining the combined effects of the Zeeman interaction 
and hybridization on the magnetocapacitance of a TI thin film device. 

Results 

We consider Dirac fermions on the surface of a 3D TI thin film which we take to be in the xy-plane. There is an 
external magnetic field perpendicular to the thin film and we take into account hybridization between the upper 
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and lower surfaces. The two-dimensional Hamiltonian for Dirac 
fermions in a magnetic field with hybridization is given by 24 ' 27,28 

7 y n x ) + (t z A z + A h )(T z . 



H — v(a r n„ 



(1 



Here, a x and a y are Pauli matrices that operate on a combined space 
of spin and surface degrees of freedom, x z = +/— denotes symmet- 
ric/antisymmetric combinations of the upper and lower surfaces of 
the TI, v denotes the Fermi velocity of the Dirac fermions, n = p + 
eA/c is the two-dimensional canonical momentum with vector 
potential A, and c is the speed of light. Moreover, we have the 

Zeeman energy A z = -gfigB, the effective Lande factor g, the Bohr 

magneton ji B , and the hybridization matrix element Af,, which 
reflects the hybridization between the upper and lower surfaces of 
the TI. We employ the Landau gauge and express the vector potential 
as A = (0, Bx, 0). The LL energies then are given by 



E? = -(A z + T z A ft ) 



E T n z A = lj2nh 2 m 2 + (A z + z z A h ) \ n^O 



(2) 



with X = + / — for electron/hole bands. Moreover, oj = v-sj eB/fi is the 
cyclotron frequency of the Dirac fermions and n is an integer. The 
main point to note about the energy spectrum, which will be crucial 
to the following discussion, is: The splitting of the energy levels 
depends on both the Zeeman interaction and the hybridization 
between the upper and lower surfaces. If either of these two mechan- 
isms is absent, no splitting will occur. Further, the energy spectrum is 
electron-hole symmetric in the absence of the Zeeman interaction, 
with a hybridization gap. We note that angular resolved photoelec- 
tron spectroscopy finds no electron-hole symmetry even in the 
absence of both Zeeman interaction and hybridization, reflecting a 
quadratic term in the Hamiltonian 39,40 . However, in our surface 
Hamiltonian in Eq. (2) no quadratic term is present as it can be 
neglected when the system is doped such that the Dirac point comes 
close to the charge neutrality point, which is the focus of our work. In 
the presence of both Zeeman interaction and hybridization, the n = 
0 LL explicitly shows a quantum phase transition. This will be 
important in the discussion of the DOS and the capacitance results. 

We consider a gated TI device in which the capacitor is formed 
between the gate and the TI thin film. The magnetocapacitance of the 
device represents the charge response in the channel as the channel 
potential is varied. In conventional devices, the magnetocapacitance 
Cq is large and can be ignored. On the contrary, in low dimensional 
devices, such as 2D electron gases, graphene, and future TI devices, it 
is the dominant capacitive contribution and therefore an important 
quantity in the design of nanoscale devices. The central expression 
for the quantum capacitance is 36 



C Q = 



edQ 



e 2 dn e 
8e F 



= e 2 D T {B), 



(3) 



where n e is the carrier concentration and t: P the Fermi energy. The 
temperature dependent DOS at a finite magnetic field, Dj{B), is 
obtained from the relation 



Dt(B) = 



dn e 



8f(e-e F ) 
de F 



D(e), 



(4) 



with the Fermi Dirac distribution/. In the limit of zero temperature 
we have D T (B) = D(t: F ) and 



(5) 



Assuming a Gaussian broadening of the LLs, the DOS per unit area is 
given by 



D(e F ) ■■ 
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where F is the width of the Gaussian distribution (zero shift). For n 
0 we have 

1 1 {e, F -El^ 2 

D - c(£f) = 2rf^rV2i exp 



2r 2 
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where = — ( A z + t z A/,) . Combining Eqs. (6) and (7), we obtain 
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At the charge neutrality point (gp = 0) the result is 



D(e c = 0) = ( -'-V— *=exp 



(A z + t z A h 



2T l 



1+2^ exp 

n = i 



which can be written as 

D( EP = 0): • 1 
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where % = ffa) 2 /2T 2 . Equation (10) reduces to the result reported in 
Ref. 36 in the absence of Zeeman interaction and hybridization. 

The magnetocapacitance Cq is plotted in Figs. 1 and 2 as a function 
of the Fermi energy (i.e., of the gate voltage). The following para- 
meters are employed 22,30,41 " 43 : g = 60, B = 3 T, A z = 5 meV, and Aj, = 
3 meV. To obtain analytical results, we choose a constant level width 
of r = 0.3 meV. We are interested in changes of the character at n P = 
0 on variation of the Zeeman interaction and hybridization relative to 
each other. It must be noted that the n = 0 LL here plays the most 
important role. Figure 1 shows that Cq is zero at £ F = 0. This occurs 
because the n = 0 LL splits only when the hybridization does not 
vanish. The n = 0 LL splits into one electron and one hole level, 
which reflects a metal to insulator transition caused by the 
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Figure 1 | Quantum capacitance as a function of the Fermi energy at 
T = 0 K, B = 3 T, and Zeeman energy 0 meV. The hybridization energy is 
0 meV (dotted lines) and 4 meV (solid lines), respectively. 
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Figure 2 | Quantum capacitance as a function of the Fermi energy at 
T = 0 K, B = 3 T, Zeeman energy 5 meV, and hybridization energy 
3 meV. 

hybridization (for A z = 0). The states are electron-hole symmetric at 
this stage. Cq can be tuned from a minimum to a maximum when the 
Zeeman energy is increased by changing the external magnetic field. 
For A z > A h both n = 0 sublevels are located in the hole region, which 
breaks the electron-hole symmetry, see Fig. 2, and represents the 
trivial to non-trivial topological insulator phase transition. To 
observe the splitting, the broadening of the LLs must be less than 
the hybridization energy. We note that not only the n = 0 LL but all 
LLs are split into two sublevels. 

In Fig. 3 we address the SdH oscillations in Cq. The following 
parameters are used: T = 0.3 meV, T = 0 K, n e = 4 X 10 15 m~ 2 , v 
= 3 X 10 5 ms _1 , and A/, = 4 meV. For low and high magnetic field 
we observe a beating pattern and a splitting of the SdH oscillations, 
respectively. Both occurs due to interference of the SdH oscillations 
at two different frequencies and therefore is a consequence of the 
splitting of the LLs for finite Zeeman interaction and hybridization. 
The beating pattern vanishes once the Zeeman interaction dominates 
the hybridization. For the chosen parameters this occurs at a mag- 
netic field of about 1.5 T. Above this value we find a well resolved 
splitting of the SdH oscillations in Fig. 3. A similar beating pattern 
and splitting of SdH oscillations in the longitudinal resistivity is seen 
in Ref. 44, where it occurs because of splitting of the LLs as a result of 
a broken inversion symmetry due to substrate effects. The sample 
investigated in this reference is thick and hybridization effects are 
ruled out. 




4 6 
Magnetic Field (Tesla) 

Figure 3 | Quantum capacitance as a function of the magnetic field at 
T = 0 K and hybridization energy 4 meV. 
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In order to understand the origin of the beating pattern and split- 
ting of the SdH oscillations quantatively, we employ the Poisson 
summation formula 



-F(0)+$>(fc)= F(n)dn + 2^2(-l) 



F(n) cos[2nkn]dn (11) 



to arrive at the following expression for the DOS at zero temperature 
which incorporates the contributions of the n s 1 LLs: 

1 h 2 0J 2 



D(e p )=Do \ 1+2 (-l)*exp 



^(4-(A z + r z A fc ) 2 ) 
n co 2 



(12) 



Here, D 0 = \ep\lnffv 2 is the zero magnetic field DOS. For T^hco, it is 
sufficient to retain only the first order term (k= 1), since the higher 
order terms are highly damped. From Eq. (12) we can confirm that 
the beating pattern is due to the interference of waves with two 
different frequencies. The oscillatory part of the k= 1 term in the 
sum for %.= +!— is 



2jc 



tfco 2 



•(A z 



2k 

tfco 2 



-(A z 



and can be expressed as 

271 

2 cos 
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This represents a wave of higher frequency whose amplitude oscil- 
lates at a lower frequency, giving rise to the beating pattern, because 
A z Ah«e, F . The amplitude of the SdH oscillations is modulated by 
4tc 



tfco 2 



;A z A h 



, , 4A Z A,, 
and nodes occur at — — — = + 0.5, : 



tfco 2 



: 1.5, where 



occurs when both 



the modulating cosine vanishes. We also note that the amplitude 

4tc 

modulation given by the term cos — ; — ^A z Ah 

[hi co 

Zeeman interaction and hybridization are present in the system. 
Moreover, the threshold magnetic field where beating is seen 
depends on the strenghts of the Zeeman interaction and hybridiza- 
tion. For magnetic fields down to about 1.5 T the beating persists. 
Above this value it is quenched to show splitting of the SdH oscilla- 
tions. As mentioned before this finding is in agreement with the 
experimental work in Ref. 44. 

Discussion 

Topological systems are typically doped away from the charge neut- 
rality point. In addition, the bulk usually dominates the surface con- 
tributions. However, thin films of topological insulators make it 
possible to access the surface contributions only 22,41 43 , typically when 
the thickness is below 5 nm. In this context, we compare our results 
with the experimental observations for the quantum magnetocapa- 
citance 45 . The latter authors study films that are thick enough that 
hybridization between the top and bottom surfaces can be ignored 
(10 nm Bi 2 Se 3 samples). To compare their results for a tilted mag- 
netic field, we next consider this situation and explore the crossover 
from external perpendicular to parallel (in-plane) magnetic field on 
the quantum capacitance. For a tilt angle 0° £ 0 £ 90° with the 
normal of the film, the magnetic field is given by (B sin 0, 0, B cos 
0). To describe the effects of the tilting we start from the two-dimen- 
sional Hamiltonian 



H — m z ya x Ky — Oy% 



(15) 



where A = (0, xB cos 0 — zB sin 0, 0) and z is the thickness of the film. 
The LL energies are obtained as 
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Figure 4 | Quantum capacitance as a function of the Fermi energy for 
different tilt angles of the magnetic field: 0° (dashed lines), 40° (red solid 
lines), 70° (blue solid lines), and 80° (black solid lines). We use B = 1 T 

and T = 25 K. 

E n ^ = ZV2nhio t , (16) 

with oj t = V\J eB cos 0/h Note the dependence of the LLs on the tilt 
angle 0 as compared to Eq. (2). 

Following the same procedure as before, we employ the energy 
spectrum to compute the quantum capacitance, which is plotted in 
Fig. 4 as a function of the Fermi energy (gate voltage) for different tilt 
angles. We focus on the broadening of the LLs and the corresponding 
effects on the SdH oscillations. We fix B = 7 T. The SdH oscillations 
are suppressed as we increase the tilt angle, since we are decreasing 
the perpendicular component of the magnetic field. For 0 — > 90°, 
when the magnetic field is almost completely aligned with the film, 
the SdH oscillations are washed out. The broadening of the LLs 
generally depends on the magnetic field strength, LL index, and 
scattering parameters. This requires a self consistent calculation, 
which usually is performed numerically. In order to carry out a 
tractable analytical calculation we choose a constant level width of 
r = 25 K. The results in Fig. 4 are consistent with the experimental 
Fig. 4 in Ref. 45. 

An experimental investigation of transport arising from the sur- 
face states in TIs is hindered by the bulk contributions due to nat- 
urally occurring defects and residual carrier doping 46 4S . In this 
regard, the low mobility of the surface states is a serious problem. 
However, recent work on TIs has shown a way to overcome this 
problem by studying transport in ultra thin samples, where a 
12 times higher surface mobility as compared to the bulk can be 
achieved 49 . It has been demonstrated that the large surface to volume 
ratio of Bi 2 Se 3 thin films effectively suppresses bulk effects and 
enables probing of the electronic transport of the surface states 50 . 
There has been a rapid improvement in the sample quality coming 
along with a reduction of film thicknesses 51 . Relevant to the present 
study, the effect of the Zeeman exchange field in TIs is significant 41,51 
with a magnitude around 50 meV due to the large effective g-factor. 
Also, recent experiments have found a hybridization gap in the ultra- 
thin regime (less than 6 nm) 22,26 as well as quantum oscillations 43 . 
Therefore, it can be expected that the effects discussed before can be 
realized in quantum capacitance experiments in ultrathin films. 

In conclusion, we have investigated the effects of Zeeman inter- 
action and hybridization on the magnetocapacitance of a TI thin film 
in a perpendicular magnetic field. To this aim, we have obtained an 
analytic expression for the DOS incorporating both the Zeeman and 
hybridization contributions. It turns out that a combination of 
Zeeman interaction and hybridization breaks the electron-hole sym- 
metry, resulting in significant effects on the magnetocapacitance at e p 



= 0. In addition, we have shown that in the presence of both Zeeman 
interaction and hybridization, the magnetocapacitance exhibits beat- 
ing at low and level splitting at high perpendicular magnetic field. 
This behaviour is explained by interference of the SdH oscillations at 
the two frequencies induced by the splitting of the LLs. Let us com- 
ment on the relevance of our findings for experiments. The cyclotron 
energy for B = 1 T is fyx>= 36 meV. Observation of the splitting of 
the LLs and the discussed consequences requires that the temper- 
ature is low and the system clean enough that disorder effects do not 
wash out the splitting. Temperature and disorder broadening thus 
must not reach the Zeeman or hybridization energies. For a tilted 
magnetic field with 0 — » 90° we find that the SdH oscillations are 
washed out. 
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